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Abstract. We consider a scenario in which primordial scalar perturbations are generated when 
complex conformal scalar field rolls down its negative quartic potential. Initially, these are the 



perturbations of the phase of this field; they are converted into the adiabatic perturbations at a 
later stage. A potentially dangerous feature of this scenario is the existence of perturbations in the 



radial field direction, which have red power spectrum. We show, however, that to the linear order 



J> ' in the small parameter — the quartic self-coupling — the infrared effects are completely harmless, 

~l_ ' as they can be absorbed into field redefinition. We then evaluate the statistical anisotropy inherent 

in the model due to the existence of the long-ranged radial perturbations. To the linear order in 
the quartic self-coupling the statistical anisotropy is free of the infrared effects. The latter show 
up at the quadratic order in the self-coupling and result in the mild (logarithmic) enhancement of 
the corresponding contribution to the statistical anisotropy. The resulting statistical anisotropy is 
a combination of a larger term which, however, decays as momentum increases, and a smaller term 
which is independent of momentum. 
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1 Introduction and summary 

The two basic properties of primordial scalar perturbations in the Universe are approximate Gaus- 
sianity and approximate flatness of the power spectrum [1]. The first property strongly suggests 
that these perturbations originate from amplified vacuum fluctuations of nearly linear (i.e., weakly 
coupled) quantum field(s): free quantum field in its vacuum state obeys the Wick theorem, the 
defining property of Gaussian random field, while linear evolution in classical backgrounds does not 
induce the non-Gaussianity. The flatness of the power spectrum calls for some symmetry behind 
it. The best known candidate is the symmetry of the de Sitter space-time under spatial dilatations 
supplemented by time translations. This is the approximate symmetry of the inflating Universe [2], 
and, indeed, the inflationary mechanism of the generation of scalar perturbations [3] produces al- 
most flat power spectrum. Another example is the symmetry (of the field equation) of the scalar 
theory with negative exponential scalar potential in flat space-time. This is the symmetry under 
space-time dilatations suppelemented by the shifts of the field. This symmetry remains approx- 
imately valid in slowly evolving, e.g., ekpyrotic [4] or "starting" [5] Universe, and the resulting 
perturbation spectrum is again almost flat 1 [6] (see also ref. [7]). 

Yet another symmetry that may be responsible for the approximate flatness of the scalar 
spectrum has been proposed in ref. [11]. This is a combination of conformal invariance and a global 
symmetry. The simplest model involves complex scalar field eft, which is conformally coupled to 
gravity and for long enough time evolves in negative quartic potential 

V(cf>) = -h 2 \cj>\ A • (1-1) 



1 There are other mechanisms capable of producing flat or almost flat scalar spectrum [8-10]. In some cases, there 
is no obvious symmetry that guarantees the flatness, i.e., the scalar spectrum is flat accidentally. 



A necessary condition for the absence of strong non-linearities at the classical level and strong 
coupling at the quantum level is 

h< 1 . 

As discussed in section 2.2, there may or may not be stronger bounds on h. The global symmetry 
in this simplest case is U(l) acting as <p — > e za 4>. One assumes that the background space-time 
is homogeneous, isotropic and spatially flat, ds 2 = a 2 (?])(dr] 2 — dx 2 ). Then, due to conformal 
invariance, the dynamics of the field 

is independent of the evolution of the scale factor and in terms of the conformal coordinates (77, x) 
proceeds in the same way as in Minkowski space-time. One begins with the homogeneous back- 
ground field Xc(v) that rolls down the negative quartic potential. Its late-time behavior is completely 
determined by conformal invariance. As we review in section 2, this rolling field produces effective 
"horizon" for perturbations 5x'- at early times the linear perturbations oscillate in conformal time 
as modes of free massless (quantum) scalar field, 2 while at late times the oscillatory behavior ceases 
to hold. The perturbations of the phase = Arg <f> freeze out at the time of the "horizon" exit, 
and their power spectrum is flat after that, 

V^=^. (1.2) 

As discussed in ref. [11], this property is a consequence of conformal and global symmetries. 

The scenario proceeds with the assumption that the scalar potential V(4>) has, in fact, a 
minimum at some large value of \cp\, and that the modulus of the field </> eventually gets relaxed to 
the minimum, see figure 1. The simplest option concerning further evolution of the perturbations 
59 is that they are superhorizon in the conventional sense by the time the conformal rolling stage 
ends. Then they remain frozen out, 3 and their power spectrum remains flat. At some much 
later cosmological epoch, the perturbations of the phase are converted into the adiabatic scalar 
perturbations; we discuss possible mechanisms responsible for that in section 2.2. These mechansims 
do not distort the power spectrum, so the resulting adiabatic perturbations have flat primordial 
power spectrum. If conformal invariance is not exact at the rolling stage, the scalar power spectrum 
has small tilt, which depends on both the strength of the violation of conformal invariance and the 
evolution of the scale factor at the rolling stage [14]. 

A peculiar, and potentially dangerous property of the model is that the modulus of the rolling 
field also acquires perturbations. Super-" horizon" modes of the modulus (i.e., radial direction) 
have red power spectrum (see section 2.3 for details), 



TWJOoc*" 1 . (1.3) 



2 An assumption here is that the rolling stage is long enough in conformal time, so that the modes of interest 
are indeed sub-" horizon" at early times. This assumption is non-trivial: the conformal rolling should occur at 
a cosmological epoch, preceeding the hot Big Bang stage, when the standard horizon problem is solved, at least 
formally, see ref. [11] and section 2. 

3 For contracting Universe, this property of superhorizon modes holds if the dominating matter has stiff equation 
of state, w > 1. This appears to be necessary for the viability of the bounce scenario anyway, see the discussion in 
refs. [12, 13]. 
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Figure 1. The scalar potential. Bullets show the evolution of the scalar field. Arrows at the end point at 
the bottom of the potential indicate that there are perturbations of the phase. 



One consequence is that there are perturbations of the energy density with red spectrum right after 
the conformal rolling stage, but before the modulus freezes out at the minimum of V{<j>). These 
are not dangerous, provided that the energy density of the field (j> is small compared to the total 
energy density at all times preceding the time when the modulus of this field settles down to the 
minimum of its potential, i.e., the cosmological evolution is governed by some other matter at that 
early epoch. In this paper we assume that this is indeed the case. 

The second consequence is that the infrared radial modes interact with the perturbations of 
the phase, and in principle may have strong effect on the latter. This is precisely the issue we 
address in this paper. To this end, we study the perturbations of the phase in the presence of 
long-ranged perturbations of the modulus and make use of the spatial gradient expansion of the 
latter. We interpret the effects emerging at the zeroth and first orders in the gradient expansion as 
the local time shift and local Lorentz boost of the background field Xc the latter property meaning 
that the rolling background is approximately homogeneous in a local reference frame different 
from the cosmic frame where the metric has the standard Friedmann-Robertson- Walker form. Our 
interpretation makes it straightforward to obtain the expression for the perturbation of the phase, 
valid to these two orders of the gradient expansion. We then show that to the linear order in h, 
the infrared effects cancel out: the perturbations of the phase remain Gaussian random and have 
fiat spectrum (1.2). Since the red power spectrum is characteristic of S\(j>\ and di(5\4>\), but not of 



higher spatial derivatives of S\(f>\, our analysis is sufficient to show that, in fact, there is no gross 
modification of the results of the linear analysis due to the effect of the infrared modes. 

The large wavelength modes of 5\cp\ are not entirely negligible, however. The modes whose 
present wavelengths exceed the present Hubble size Hq induce statistical anisotropy in the pertur- 
bations of the phase 56, and hence in the resulting adiabatic perturbations. To the linear order in 
h, the statistical anisotropy is generated at the second order in the gradient expansion, and hence 
it is free of the infrared effects. The infrared modes show up at the second order in h, and lead to 
mild enhancement of the statistical anisotropy at this order. Accordingly, the power spectrum of 
the adiabatic perturbation £ has the following form, 

P c (k) = V (k) fl + d-h-^-- kikj*>ij ~C2-h 2 - (feu) 2 ) . (1.4) 

In the first non-trivial term, Wij is a traceless symmetric tensor of a general form with unit nor- 
malization, WijWij = 1, k is a unit vector, k = k/fe, and c\ is a constant of order 1 whose actual 
value is undetermined because of the cosmic variance. As we alluded to above, the deep infrared 
modes are irrelevant for this term. In the last term, u is some unit vector independent of Wij, and 
the positive parameter C2 is logarithmically enhanced due to the infrared effects. This is the first 
place where the deep infrared modes show up. Clearly, their effect is subdominant for small h. 

The statistical anisotropy encoded in the second term in (1.4) is similar to that commonly dis- 
cussed in inflationary context [15], and, indeed, generated in some concrete inflationary models [16]: 
it does not decay as momentum increases and has special tensorial form (ku) 2 with constant u. On 
the other hand, the first non-trivial term in (1.4) has the general tensorial structure and decreases 
with momentum. The latter property is somewhat similar to the situation that occurs in cosmolog- 
ical models with the anisotropic expansion before inflation [17]. Overall, the statistical anisotropy 
(1.4) may be quite substantial, since there are no strong bounds on h at least for one mechanism 
of conversion of the phase perturbations into adiabatic ones, as discussed in section 2.2. 

Interestingly, the overall statistical anisotropy is a combination of a larger (order 0(h)) term 
which, however, decays as k increases, and a smaller (order 0(/i 2 logA), where A is the infrared 
cutoff) term which is independent of k. We consider this feature as a potential smoking-gun 
property of our scenario. 

It is worth noting that the non-linearity of the scalar potential gives rise to the non-Gaussianity 
of the perturbations of the phase 56, and hence the adiabatic perturbations in our scenario, over 
and beyond the non-Gaussianity that may be generated at the time when the phase perturbations 
get reprocessed into the adiabatic perturbations. In view of the result outlined above, this non- 
Gaussianity is not plagued by the infrared effects at the first non-trivial order in h, and must be 
fully calculable at this order. We do not consider the non-Gaussianity in this paper, since the 
derivative expansion approach we employ here is useless in this regard. 

The paper is organized as follows. In section 2 we review the linear analysis of the model. 
In section 3 we study the effect of infrared modes of the modulus 5\<f>\ on the perturbations of the 
phase 56 at the leading and subleading orders of the gradient expansion, and to the linear order 
in h. Statistical anisotropy, which is generated at the next order, is analysed in section 4. We 
conclude in section 5. 



2 Linear analysis 

At the conformal rolling stage, the dynamics of the scalar field is governed by the action 



where the scalar potential V(4>) is negative and has conformally invariant form (1.1). In terms of 
the field x = a< />) the action in conformal coordinates is the same as in Minkowski space-time, 

S[ X ] = j d 3 xd v [v^d^d uX + h 2 \x\ 4 ] ■ 
The field equation is 

rrd ll a v x-2h 2 \x\ 2 x = o. (2.1) 

Spatially homogeneous solutions approach the late-time attractor 

*w = 55^ • (2 ' 2) 

where ??* is an arbitrary real parameter, and we consider real solution, without loss of generality. 
We take the solutuion (2.2) as the background. The meaning of the parameter 77* is that the field Xc 
would run away to infinity as r/ — > 77* , if the scalar potential remained negative quartic at arbitrarily 
large fields. 

2.1 Perturbations of phase 

At the linearized level, the perturbations of the phase and modulus of the field eft decouple from each 

other. Let us begin with the perturbations of the phase, or, for real background (2.2), perturbations 

of the imaginary part 

Im x 

They obey the linearized equation, 

(6x2)" ~ dA 8x2 ~ 2h\ 2 c 5x2 = 0, (2.3) 

where prime denotes the derivative with respect to r\. Explicitly, 

(S X 2)" - dA 5 X 2 - ( 2 _ )2 5x2 = 0. (2.4) 

Let k be conformal momentum of perturbation. An important assumption of the entire scenario is 
that the rolling stage begins early enough, so that there is time at which the following inequality 
holds: 

fc(r?*-r/)»l- (2.5) 

Since the momenta k of cosmological significance are as small as the present Hubble parameter, 
this inequality means that the duration of the rolling stage in conformal time is longer than the 



conformal time elapsed from, say, the beginning of the hot Big Bang expansion to the present epoch. 
This is only possible if the hot Big Bang stage was preceded by some other epoch, at which the 
standard horizon problem is solved; the mechanism we discuss in this paper is meant to operate at 
that epoch. We note in passing that the latter property is inherent in most, if not all, mechanisms 
of the generation of cosmological perturbations. 

Equation (2.4) is exactly the same as equation for minimally coupled massless scalar field in 
the de Sitter background. Nevertheless, let us briefly discuss its solutions. At early times, when 
the inequality (2.5) is satisfied, 6x2 is free massless quantum field, 



/d 3 k ( (-\ 
; — 1= [5x1 '( k, x, ri)A\ + h.c. 



d 3 k 

whose modes are (we keep the dependence on x for future convenience) 

^" ) (k,x,r ? )=e ikx -^. (2.6) 

Here A^ and A^ are aniihilation and creation operators obeying the standard commutational re- 
lation, [.Aio-Ajj/] = <5(k — k'). It is natural to assume that the field 5x2 is initially in its vacuum 
state. 

The rolling background Xc(jl) produces an effective "horizon" for the perturbations 6x2- The 
oscillations (2.6) terminate when the mode exits the "horizon", i.e., at k(j)* — rj) ~ 1. The solution 
to eq. (2.4) with the initial condition (2.6) is 

5x { 2 - ] (k, x, n) = e* kx -^* • F(k, n*-n), (2.7) 

where 

/^7T (i) 



F{KZ) = -^-kiHl)> 2 {ki) (2. 

i o +ri£i I— Ion \ra I Tmnr»+i(-ivi It 

3/2 



and H^ /2 is the Hankel function. In the late-time super-" horizon" regime, when k(rj* — rj) <C 1, one 
has 



Hence, the super-" horizon" perturbations of the phase 56 = 5x2 /Xc are time-independent, 

_ fo2 (x, 7?) _ f d 3 k ikx _ ifc „. j , (0 im 

() ~~XM~~ J 4^^ e A * + h - C - (2 ' 10) 

This expression describes Gaussian random field (cf. ref. [18]) whose power spectrum is given by 
(1.2). 

2.2 Converting perturbations of phase into adiabatic perturbations 

For the sake of completenetss, let us briefly discuss two possible ways to reprocess the perturbations 
of the phase 59 into the adiabatic perturbations. 

One way is to make use of the curvaton mechanism [19] specified to pseudo-Nambu-Goldstone 
curvaton [20]. Namely, let us assume that 8 is actually pseudo-Nambu-Goldstone field. Generically, 
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Figure 2. The scalar potential in the pseudo-Nambu-Goldstone scenario. 

conformal rolling ends up at a slope of its potential, see Fig. 2. The field 9 together with its 
perturbations remains frozen out until the time at the hot Big Bang epoch when the Hubble 
parameter becomes of the order of the mass parameter of this field. Then 9 starts to oscillate 
about the minimum of its potential. Provided this field interacts with conventional particles, its 
oscillations eventually decay, and their energy density is transfered to the cosmic plasma. At 
that epoch, the perturbations of the energy density of the field 8 are converted into the adiabatic 
perturbations. The shape of the resulting adiabatic power spectrum is the same as that of the 
initial spectrum of 59, i.e., it is flat in the case of exact conformal invariance at the rolling stage. 
Generically, the amplitude of the adiabatic perturbations is of order 



h 



~ r- 



2ir9 e 



where r is the ratio of the energy density of the field 9 to the total energy densiy at the time the 
^-oscillations decay, and 9 C is the distance from the landing point of rolling to the minimum of the 
potential. 

The pseudo-Nambu-Goldstone mechanism produces the non-Gaussianity of local form in the 
adiabatic perturbations. For generic values of the phase at landing, 9 C ~ vr/2, non-observation of 



the non-Gaussianity [1] implies 

r > 10" 2 , 

so that the correct scalar amplitude is obtained for 

h<10~ 2 . (2.11) 

Thus, the scenario with the pseudo-Nambu-Goldstone mechanism of conversion of the phase per- 
turbations into the adiabatic ones is viable for quite small scalar self-coupling only. 

Another possibility is the modulated decay of heavy particles [21-23]. One assumes that the 
phase field 9 interacts with some heavy particles in such a way that the masses and/or widths of 
the latter depend on 9, 

M = M + e M 9 and/or V = T Q + e T 9 . (2.12) 

One assumes further that these particles survive at the hot Big Bang epoch until they are non- 
relativistic and dominate the cosmological expansion. When these particles decay, the perturbations 
in 9, and hence in M and/or V, induce adiabatic perturbations, 

> 5M e M 59 , . > 5T e r 59 

C ~ -irr = tt 7T and/or C ~ — - = — — . 

M M + € M 9 C ' r r + e r # c 

The shape of the adiabatic power spectrum is again the same as that of the initial power spectrum 
of 59. 

The modulated decay mechanism also induces non-Gaussianity in the adiabatic perturbations. 
However, once the dependence of the mass/width on 9 is linear, as written in (2.12), the induced 
non-Gaussian part of the adiabatic perturbations is of order 

(SM/M) 2 , (5T/T) 2 ~ C 2 • 

In other words, irrespectively of the value of the coupling h, the non-Gaussianity parameter is fairly 
small (see refs. [22, 23] for details), 

Inl ~ 1 , 

in comfortable agreement with the existing limit [1]. Thus, the modulated decay mechanism by 
itself does not imply any bound on h. 

2.3 Perturbations of modulus 

Let us now come back to the conformal rolling stage and consider the radial perturbations — pertur- 
bations of the modulus of the field X-> or > with our convention of real backgronud Xc ; perturbations 
of the real part 

Re x 

At the linearized level, they obey the following equation, 

6 
(5xi)" - didi 5xi - Qx 2 c 5xi = (5xi)" ~ didi 5xi - -, Tn^Xi = . 

iv* - vr 



Its solution that tends to properly normalized mode of free quantum field as £(77* — 77) — )• 00 is 



5X1 = e ikx ~*"* • t^V^V^^S IHrt* ~ V)] ■ Ac + fc.c , 



4vr V 2 5 /2 

where -Bk; B^ is another set of annihilation and creation operators. At late times, when £(77* — 77) <C 
1 (super-" horizon" regime), one has 

^i=e^-^- i j 7 , fe5/2( - 1 _ 7?)2 -i k + /,.c. 

Hence, the super-" horizon" perturbations of the modulus have red power spectrum (1.3). 

The dependnce 6x1 ^ (v* ~ r ?)~ 2 i s naturally interpreted in terms of the local shift of the "end 
time" parameter 77*. Indeed, with the background field given by (2.2), the sum Xc + Sxi, i- e -> the 
radial field including perturbations, is the linearized form of 

Xc[r?*(x) - 77] = 7 , (2.13) 

%*(x) -77] 

where 

??*(x) = ??* + <5t7*(x) (2.14) 

and 

So, the infrared radial modes modify the effective background by transforming the "end time" 
parameter 77* into random field that slowly varies in space, 4 as given in eqs. (2.13), (2.14). Clearly, 
this observation is valid beyond the linear approximation: once the spatial scale of variation of 
Xi(x, 77) exceeds the "horizon" size, spatial gradients in eq. (2.1) are negligible, and the late-time 
solutions to the full non-linear field equation have locally one and the same form (2.2), modulo 
slow variation of 77* in space. 

It is worth noting that the infrared modes contribute both to the field 5rj*(x) itself and to its 
spatial derivative. The contribution of the modes which are superhorizon today, i.e., have momenta 
k < Ho, to the variance of the latter is given by 

(d i ?7*(x)d j ?7*(x)) fe <Ho = <V 4^- I -J- = °ij ' 4^ log X ' ^ ' 

where A is the infrared cutoff which parametrizes our ignorance of the dynamics at the beginning 
of the conformal rolling stage. 

3 Effect of infrared radial modes on perturbations of phase: first order in h 

The main purpose of this paper is to understand how the interaction with the infrared radial modes 
affects the properties of the perturbations of the phase 56. To this end, we consider perturbations of 



4 There are corrections to eq. (2.13) of order didjr]*(x.) and [dirj t (x)] 2 , see sections 4.1 and 4.2. 



the imaginary part 5x2 , whose wavelengths are much smaller than the scale of the spatial variation 
of the modulus. Because of the separation of scales, perturbations 5x2 can still be treated in the 
linear approximation, but now in the background (2.13). 

Since our concern is the infrared part of 7/*(x), we make use of the spatial gradient expansion, 
consider, for the time being, a region near the origin and write 

r/*(x) = 77* (0) - ViXi + . . . , (3.1) 

where 

Vi = -<9j?7*(x)| x= o , 

and dots denote higher order terms in x. Importantly, the field didjrj*(x) has blue power spectrum, 
unlike ry*(x) and c^r/^x), so the major effect of the infrared modes is accounted for by considering 
the two terms of the gradient expansion written explicitly in (3.1). In this section we work at this, 
first order of the gradient expansion. Furthermore, we assume in what follows that 

|v| < 1 , (3.2) 

and in this section we neglect corrections of order v 2 . The expansion in |v| is legitimate, since 
the field v(x) has flat power spectrum, so the fluctuation of v is of order h 2 \ log A|, where A is the 
infrared cutoff, and it is small for small h and not too large | log A|. In other words, the expansion 
in |v| is the expansion in h, modulo infrared logarithms. We postpone to section 4.2 the analysis 
of the leading effect that occurs at the order v 2 . 

Keeping the two terms in (3.1) only, we have, instead of eq. (2.4), 

{5X2)" - didi 5 X 2 - r , m 2 t^5x2 = . (3.3) 

[77* (0) - vx- i]\ 2 

We observe that the denominator in the expression for the background field 

h[r]*(0) — rj — vxj 

contains the combination 77* (0) — (77 + vx). We interpret this as the local time shift and Lorentz 
boost of the original background (2.2): the effective background is homogeneous in a reference 
frame (in conformal coordinates) other than the cosmic frame where the metric is spatially homo- 
geneous. Note that the field (3.4) is a solution to the field equation (2.1) in our approximation. 
Our interpretation suggests that the solutions to eq. (3.3) can be obtained by time translation and 
Lorentz boost of the original solution (2.7), (2.8). Indeed, it is straightforward to see that to the 
first order in v, the solution to eq. (3.3) obeying the initial condition (2.6) is 

SxP (k, x, 77) = e il(x+w)-i<iv* (0) . F ( ?j ^ (0) - 77 - vx) , (3.5) 

where the function F is still defined by (2.8), the Lorentz-boosted momentum is 

q = k + kv , q = |q| = k + kv , (3.6) 
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and it is understood that terms of order v 2 must be neglected. The solution (3.5) can be cast into 
the following form, 

8 X i-\Kx,v) = e ikx -^*( x )- kv M*)-.] . F(g,ifc(x) - r,) , (3.7) 

where 

Vi = -9j77*(x) . 

This form, still valid to the first order in the gradient expansion and to the linear order in v, does 
not make any reference to the origin of the coordinate frame and can be used for arbitrary x. We 
consider corrections of order didjT]*(x.) and v 2 to this solution in sections 4.1 and 4.2, respectively. 
We find from eqs. (3.4) and (3.5) (or, equivalently, (2.13) and (3.7)) that the perturbations of 
the phase again freeze out as /c[r/*(x) — rj] — > 0, now at 



Xc(x) J y/k A-K 6 l 2 q 



(3.8) 



where the relation between k and q is still given by eq. (3.6). We recall that 17* (x) here is a 
realization of infrared random field. The formula (3.8) implies that to the first order of the gradient 
expansion we limit ourselves in this section, the properties of the random field 59 are the same as 
those of the linear field (2.10). Consider, e.g., the two-point correlation function 

/A'ih, u2 
-, ^—-e ik(xi - X2) - ifc[ "* (xi) - , '* (x2)] +c.c. (3.9) 

k 167r d g 2 

Infrared modes of 77* (x) we consider have momenta much lower than k, i.e., the spatial scale of 
their variation much exceeds |xi — X2I. Hence, to the first order in the gradient expansion we have 
77*(xi) — 77* (x 2 ) = — v(xi — X2), and v is independent of x. Therefore, the two-point function is 

(W( X1 ) W(x a )> = J — ^e^^) + c.c. 

We now change the integration varible from k to q, recall that the integration measure d s k/k is 
Lorentz-invariant, and obtain 

(W(xi) <W(x 2 )> = h 2 J J^e^-*^ + c.c. 

This is precisely the two-point correlation function of the linear field (2.10). 

The latter argument is straightforwardly generalized to multiple correlators: for a given re- 
alization of the random field 7?*(x), they are all expressed in terms of the two-point correlation 
function (3.9). In other words, the infrared effects are removed by the field redefinition, 

i q = e-*"*( )J-i k , (3.10) 

where k and q are still related by (3.6). The operators Aq, Aq obey the standard commutational 
relations, while in our approximation, the field (3.8), written in terms of these operators, coincides 
with the linear field (2.10). We conclude that the infrared radial modes are, in fact, not particularly 
dangerous, as they do not grossly affect the properties of the field 59. 
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4 Statistical anisotropy 

4.1 First order in h 

Let us continue with the analysis at the first order in h. In this approximation, the non-trivial 
effect of the large wavelength perturbations 5r]* (x) on the perturbations of the phase, and hence on 
the resulting adiabatic perturbations, occurs for the first time at the second order in the gradient 
expansion, i.e., at the order didji]*. Let us concentrate on the effect of the modes of <5r/* whose 
present wavelengths exceed the present Hubble size. We are dealing with one realization of the 
random field 5?7*, hence at the second order of the gradient expansion, didji]* is merely a tensor, 
constant throughout the visible Universe. In this section we calculate the statistical anisotropy 
associated with this tensor. 

To this end, we make use of the perturbation theory in didji]*. In the first place, we have to 
find the background, since the function (2.13) is no longer the solution to the field equation (2.1) 
at the second order in the gradient expansion. We write 

Xc = X? ) +X? ) , (4.1) 

where 

yW = - (4 2) 

(2) 

and Xc = 0{didjT] if ). We plug the expression (4.1) into eq. (2.1), linearize in didjrj* and obtain 
the following equation for the correction, 

xi 2)// -^-^xi 2) =-T7^^, (4-3) 



where we again neglect terms of order v 2 . Clearly, the correction to the background depends only 
on the scalar didirj*, so it does not yield statistical anisotropy. Nevertheless, we keep this correction 
in what follows. 

The general solution to eq. (4.3) is 

X? )= Clh .-,)3 + I -^- ? + i- 8i8j „. 

The first term in the right hand side is irrelevant in the super-" horizon" regime, the second term 
is merely a shift of 77* in the leading order expression (4.2), so the non-trivial correction is given by 
the third term. The combination entering eq. (2.3) for the perturbations of the imaginary part is 
now given by (we consistently work at the linear order in didji]*) 

2h 2 X 2 = - + - didiV * . (4.4) 

(r?*(x)-?7) 2 3^-r/ 

Note that we can set 77* = const in the second term, since we neglect corrections of order didj7]*-dkr]*- 
Let us now obtain the solution to eq. (2.3) to the first order in didjij^. The initial condition 
is still given by eq. (2.6). We search for the solution in the following form, 

6xi~\k,x,rj) = e *kx-^(xMkvMx)-^] . [p(^( x ) _ v ) + F^(q, v ^) - 77) 
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where the leading term F is again defined by (2.8) and F^ 2 > is proportional to didj?]*. We expand 
eq. (2.3) in didjrj* and obtain to the linear order 



F m " + k 2 F m _ 2 F ( 2 ) = a . a .^ . 5 + kj. g. a . r ) (45) 



where £ = 77* — 77, 



s = _*,+« + £,. ( «) 

r = -2F S - 2i^|« = - Ae«* . (4 .7) 

After calculating the right hand side of eq. (4.5) we set q = k and v = 0, since we neglect the 
terms of order didjT]* ■ dtr]*. Note that the last term in (4.6) comes from the correction to the 
background, see eq. (4.4), and that the last expression in (4.7) is obtained by using the explicit 
form of the Hankel function H\L. 

The soultion F^ 2 > should vanish as r/* — r\ — > oo. Hence, it is given in terms of the retarded 
Green's function (recall that £' > £ corresponds to r/ < r/) 

Gfo?) = ^~ ■ e(C " • [Js/2(kS)N 3/2 (k?) - N 3/2 (kOJ 3/2 (k£')] , (4.8) 

where J3/2 and ^3/2 are the Bessel functions. Namely, 

dg G(C, O [did iV * ■ 5(0 + kikj didtf* ■ T(0] • (4.9) 

We are interested in the behavior of this solution in the super-" horizon" regime, fc£ — > 0. Since 
the most singular behavior of S and T at small £ is £ -2 , the first term in square brackets in (4.8) 
gives finite contribution to the integral (4.9) as £ — > 0, while the second term yields F^ 2 '(£) oc £ _1 . 
Hence, the behavior of the correction F( 2 \£) at small £ is the same as that of the leading term F(£), 
so the phase perturbations 59, the correction included, are constant in time in the super-" horizon" 
regime. This is a consistency check of our approach. 

Since we treat didjf]* as constant over the visible Universe, the power spectrum of the phase 
perturbations contains a part linear 5 in didj?]*. This linear part comes from the interference of F^ 2 ' 
and F. In the super-" horizon" regime, F is pure imaginary, see (2.9). Hence, we are interested in 
the behavior of the imaginary part of F^ 2 > as £ — > 0. This behavior is straightforward to obtain. 
Namely, for T-term in (4.9) we write 

Im F^ T \^ = -^N 3/2 (kt;) • j™ dH'^J 3/2 {kZ>) (~^f-) • hk 3 && 



5 Had we considered the perturbations Sr/, whose wavelengths are much shorter than the present Hubble size, it 
would be appropriate to perform ensemble averaging. In that case the linear part would average out. 
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Performing similar calculation for S'-term in (4.9), we obtain that in the super-" horizon" regime 

F+fW -iS^)( 1 -a-^*** + s^)- 

The two non-trivial terms in parenthesis give the correction to the power spectrum of the phase 
perturbations due to the radial modes whose whavelengths exceed the present Hubble size. The 
same correction is characteristic of the adiabatic perturbations, so we have finally 



V C = A ( 



1 - t ■ -7/ f didjT]* - -Sijdkdk-r)* j 



(4.10) 



where the adiabatic amplitude A^ is independent of k within our aproximations. Notably, the 
angular average of the correction vanishes, so we are dealing with the statistical anisotropy proper. 
Neither the magnitude nor the exact form of the tensor didjf]* — (l/3)8ijd k dkr)* can be un- 
ambiguously predicted because of the cosmic variance. To estimate the strength of the statistical 
anisotropy, let us consider the variance 

{{ddr}*) 2 ) = ddidjT]* - -8ijd k d k r)A ■ (djdj-77* - -8ijd k d k r]A) k < Ho , 

where the notation reflects the fact that we take into account only those modes whose present 
wavelengths exceed the present Hubble size. We make use of (2.15) and obtain 

//an ^2\ 9h2 f dk 2fc4 3h2 rj2 tA ^\ 

«*>*>> = IS? J k < Ho ¥ ■ — * 1^ ' ■ (411) 

In this way we arrive at the first non-trivial term in (1.4). Higher orders in the gradient expansion 
give contributions to the statistical anisotropy which are suppressed by extra factors of Ho/k. 

4.2 Order h 2 : contribution of deep infrared modes 

Let us now turn to the statistical anisotropy at the second order in h. The major contribution 
at this order is proportional to dirj^dji]^ = V{Vj. Indeed, the overall time shift 77* (0) is irrelevant, 
while the terms involivng higher derivative combinations like didjr]*dirj* are suppressed by powers 
of Ho/k, cf. (4.10), (4.11). Furthermore, the "velocity" v is enhanced, albeit only logarithmically, 
by the deep infrared effects, see (2.16). Hence, to extract the major contribution to the statistical 
anisotropy at the order h 2 , we use the two terms of the derivative expansion, explicitly written in 
(3.1). 

Once the terms of order v 2 are not neglected, the function (3.4) is no longer a solution to the 
field equation (2.1). Instead, the solution is (keeping two terms only in the gradient expansion) 

Xc = 7 ; * T • (4.12) 

fry to* (0) -77-vxj 

So, it is appropriate to study the solutions to eq. (2.3) in this background. It is straightforward to 
see that the solution that obeys the initial condition (2.6) is 

5xi~\k,x,ri) =e i9 ll 7(x il +u " )+iqTxT - i9 ™* (0) • F [9,7(77,(0) - r/ - vx)] , (4.13) 
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where the Lorentz-boosted momentum is, as usual, q\\ = 7(^1 + kv), q T = k T , q = j(k + knv), 
7 = (1 — t; 2 ) -1 ' 2 , and notations || and T refer to components which are parallel and normal to 
v, respectively. The interpretation of this solution is again that it is the Lorentz-boosted solution 
(2.7). 

According to the scenario discussed in this paper, the phase perturbations 50 freeze out at the 
hypersurface 

7? = 77* (0) — vx = 77*(0) — vx\\ 

and then stay constant in the cosmic time rj. It follows from (4.12) and (4.13) that at this hyper- 
surface and later, the perturbations of the phase are 

<J0(x) = I ^-4r e^ll^+^^A, + h.c. , (4.14) 

where the operators Aq are defined by (3.10) and obey the standard commutational relations, and 
we have omitted irrelevant constant phase factor. We see that to the order v 2 , the effect of the 
deep infrared modes does not disappear: it is encoded in the factor 7 -1 in the first term in the 
exponent. It is clear that the momentum of perturbation labeled by q is actually equal to 

p = (7 -1 9||,q T ) • 

Accordingly, the power spectrum (omitting the correction discussed in section 4.1) is given by 

d 3 p h 2 jd 3 p 



r se (p) 



Airp 3 16vr 3 [(-y^||) 2 + (p T ) 2 ] 3 / 2 ' 



The same power spectrum characterizes the adiabatic perturbations in our scenario, and changing 
the notation p — > k we obtain finally 

p 3 ( 3 (kv) 2 

^^ = AC [( 7 fc||) 2 + (k^) 2 ]3/2 = A ({ 1 - 2-tf- 

Hence, we have arrived at the last term in (1.4). Again, neither the direction of v nor its length 
can be unambiguously calculated because of cosmic variance; recall, however, that the value of |v|, 
and hence the parameter C2 in (1.4), is logarithmically enhanced due to the infrared effects, see 
(2.16). 

In the case of the pseudo-Nambu-Goldstone mechanism of conversion of the phase perturba- 
tions into the adiabatic ones, both terms in the statistical ansotropy (1.4) are small because of the 
bound (2.11). On the other hand, the effect may be stronger in the case of the modulated decay 
mechanism. 

5 Conclusion 

We conclude by making a few remarks. 

First, our mechanism of the generation of the adiabatic perturbations can work in any cos- 
mological scenario that solves the horizon problem of the hot Big Bang theory, including inflation, 
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bouncing/cyclic scenario, pre-Big Bang, etc. In some of these scenarios (e.g., bouncing Universe), 
the assumption that the phase perturbations are superhorizon in conventional sense by the end of 
the conformal rolling stage may be non-trivial. It would be of interest to study also the opposite 
case, in which the phase evolves for some time after the end of conformal rolling. 

Second, we concentrated in this paper on the effect of infrared radial modes, and employed 
the derivative expansion. The expressions like (3.8), which we obtained in this way, must be used 
with caution, however. Bold usage of (3.8) would yield, e.g., non- vanishing equal-time commutator 
[#(x), (9(y)], which would obviously be a wrong result. The point is that the formula (3.8) is valid 
in the approximation v = const; with this understanding, the equal-time commutator vanishes, as 
it should. 

Finally, the non-linearity of the field equation gives rise to the intrinsic non-Gaussianity of 
the phase perturbations and, as a result, adiabatic perturbations. The non-Gaussianity emerges 
at the order 0(h 2 ), and may be sizeable for large enough values of the coupling h. The form of 
the non-Gaussianity is rather peculiar in our scenario. Unlike in many other cases, the three-point 
correlation function vanishes, while the four-point correlation function of 50 (and hence of adiabatic 
perturbations) involves the two-point correlator of the independent Gaussian field Sr/*. In view of 
the result of section 3, correlation functions of 50 should be infrared-finite, at least at the order 
0(h 2 ). This result shows also that the analysis of the non-Gaussianity requires methods beyond 
the gradient expansion. We leave this important issue for the future. 
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